OPEN 3 ACCESS Freely available online 



•0-PLOS I ONE 



g-Jitter Mixed Convective Slip Flow of Nanofluid past a (8) 
Permeable Stretching Sheet Embedded in a Darcian crossMaA 
Porous Media with Variable Viscosity 

Mohammed J. Uddin\ Waqar A. Khan^% Norsarahaida S. Amin^ 

1 Department of Mathematics, American International University- Bangladesh, Banani, Dhaka, Bangladesh, 2 Department of Engineering Sciences, PN Engineering College, 
National University of Sciences and Technology, Karachi, Pakistan, 3 Department of Mathematical Sciences, Faculty of Science, Universiti Teknologi Malaysia, Johor, 
Malaysia 



Abstract 

The unsteady two-dimensional laminar g-Jitter mixed convective boundary layer flow of Cu-water and Al203-water 
nanofluids past a permeable stretching sheet in a Darcian porous is studied by using an implicit finite difference numerical 
method with quasi-linearization technique. It is assumed that the plate is subjected to velocity and thermal slip boundary 
conditions. We have considered temperature dependent viscosity. The governing boundary layer equations are converted 
into non-similar equations using suitable transformations, before being solved numerically. The transport equations have 
been shown to be controlled by a number of parameters including viscosity parameter, Darcy number, nanoparticle volume 
fraction, Prandtl number, velocity slip, thermal slip, suction/injection and mixed convection parameters. The dimensionless 
velocity and temperature profiles as well as friction factor and heat transfer rates are presented graphically and discussed. It 
is found that the velocity reduces with velocity slip parameter for both nanofluids for fluid with both constant and variable 
properties. It is further found that the skin friction decreases with both Darcy number and momentum slip parameter while 
it increases with viscosity variation parameter. The surface temperature increases as the dimensionless time increases for 
both nanofluids. Nusselt numbers increase with mixed convection parameter and Darcy numbers and decreases with the 
momentum slip. Excellent agreement is found between the numerical results of the present paper with published results. 
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Introduction 

The presence of temperature gradient and gravitational field 
yields convective flows in clear as well as porous media. This type 
of flow has a significant effect on the homogenous melt growth of 
semiconductor or metal crystals on earth-bound conditions. In 
space, the gravity effect is suppressed remarkably and hence 
buoyancy effect also reduces. However, microgravity environment 
is helpful in suppressing convective flows. The effect of g-Jitter 
which originates from crew motions, mechanical vibrations 
(pumps, motors, excitations of natural frequencies of spacecraft 
structure), spacecraft maneuvers, atmospheric drag and the earth's 
gravity gradient (Li [1]) have shown to make it difficult to realize a 
diffusion-controlled growth from melts in microgravity (Lehoczky 
[2]). This adverse effect of g-jitter has been the subject of many 
investigations (see, for instance, Chen and Saghir [3]). A number 
of researchers studied g-Jitter convective flow in various aspects. 
Wasu and Rajvanshi [4] studied unsteady mixed convection flow 
under the influence of gravity modulation and magnetic field. The 
g-Jitter effects on viscous fluid flow and porous medium have also 
been investigated by Rees and Pop [5,6], Chen and Chen [7]. 

Transport phenomena associated with nanofluids flow have 
received attention of many investigators due to their diverse 
applications in many fields such as the delivery of nanodrug 



(Kleinstreuer et al. [8]), thermal management (Nkurikiyimfura et 
al. [9]), solar collectors (Said et al. [10]), transportation, the 
environment and national security (Varjavalu et al. [11] and 
Wong and Leon [12]). According to Krajnik et al. [13], nanofluids 
can be optimized during manufacture using sheet processing. 
Many superior lubricants as well as thermal working fluids may be 
developed for applications in aerospace, energy systems, medical 
engineering etc. Available mathematical models for boundary 
layer flow of nanofluids are (i) Buongiorno [14] model which 
incorporates Brownian motion and thermophoresis effects and (ii) 
the Tiwari and Das model [15] which can be used to analyze the 
behavior of nanofluids taking into account the solid volume 
fraction. Several authors have used these two models to study 
various flow phenomena past various types of geometries subject to 
various initial and boundary conditions. As an example, 
Kuznetsov and Nield [16] obtained similarity solution for natural 
convection flow of a nanofluid along a vertical plate. The effects of 
Brownian motion and thermophoresis on boundary layer flow of 
nanofluids past a flat surface in a porous medium was studied by 
Nield and Kuznetsov [17]. Chamkha et al. [18] investigated the 
natural convective flow of nanofluid past an isothermal sphere 
placed in a Darcy porous medium. Khan and Aziz [19] studied the 
double-diffusive natural convection from a vertical plate to a 
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porous medium saturated with a binary base fluid containing 
nanoparticles. Analytical and numerical analysis of heat and mass 
transfer of nanofluid was carried out by Kameswaran et al. [20] . 
Finite difference solution of transport equations due to nanono- 
fluid flow past a vertical plate taking into account double 
stratification was presented by Ibrahim and Makinde [21]. 
Cimpean and Pop [22] investigated nanofluid mixed convection 
in a porous medium channel. Nanofluid convection flows in 
porous media have received attention motivated by materials 
processing and solar energy collector applications. Gorla et al. [23] 
used a Blottner difference method to simulate mixed convection 
flow from a vertical wedge in a porous medium saturated with a 
nanofluid. Beg et al. [24] used finite difference numerical method 
to study nanofluid flow, heat and mass transfer in a porous 
medium. Further studies of nanofluid convection transport in 
porous media have been reported by Gorla and Chamkha [25] for 
non-isothermal effects, Yasin et al. [26] for heat generation effects, 
Kuznetsov [27] for non-oscillatory and oscillatory bioconvection, 
Murthy et al. [28] for magnetic effect on thermally stratified 
medium. Beg et al. [29] used a homotopy analysis method to 
simulate nanofluid free convection from a spherical body in a 
Darcian porous medium. Chamkha and Rashad [30] studied 
steady laminar natural convection boundary layer flow over a 
permeable vertical cone in a porous medium saturated with a 
nanofluid. Very recently, Uddin et al. [31], present a mathemat- 
ical model of nonlinear radiative hydromagnetic thermo-solutal 
nanofluid convection slip flow in saturated porous media. 
Relatively few studies of flow with multiple slip effect in porous 
media have been communicated. 

Much research has been focused on heat and fluid flow at 
micro-scale and nano-scales because of the applications in micro- 
electro-mechanical systems and nano-electro-mechanical systems 
with the inclusion of velocity and temperature slip boundary 
conditions at the wall, as for this type of flow convectional no slip 
conditions yield unrealistic results. In the slip region, the Navier- 
Stokes and energy equations can still be used provided the velocity 
slip and temperature slip at the walls are taken into account 
simultaneously (Karniadakis et al. [32]). Bocquet and Barrat [33] 
illustrated the mechanisms of surface slip (from nano to micro 
scales) and heat transfer on the interface. Various aspect of fluid 
flow and heat transfer of nanofluid slip flows past various 
geometries were investigated by many authors (Robert et al. 
[34], Wang [35], Zheng [36], Wu [37], Noghrehabadi et al. [38], 
Matthews and Hill [39], Das [40], Bhattacharyya et al. [41]). 

It seems that most of the published papers on nanofluid 
considered viscosity as constant. However, in reality it depends on 
temperature. Heat generated by internal friction and the corre- 
sponding rise in the temperature affects the viscosity of the fluid, so 
that the fluid viscosity no longer is considered constant (Mukho- 
padhyay and Layek [42]). Prasad et al. [43] study the effect of 
variable fluid properties on the hydromagnetic flow and heat 
transfer over a stretching sheet. The effect of temperature 
dependent viscosity and thermal conductivity on flow, heat and 
mass transfer flow was investigated by Hamad et al. [44] . Pal and 
Mondal [45] investigated the effect of variable viscosity on mixed 
convection flow and heat transfer along a stretching sheet a in a 
non-Darcy porous medium. Hence, in order to get the real picture 
of the flow, heat and mass transfer characteristics in the nanofluids, 
it is required to take into the effect of the temperature-dependent 
viscosity of the base fluid. This motivates the present study. 

In the present work, we investigate the effects of linear 
hydrodynamic slip, thermal slip and temperature dependent 
viscosity on g-Jitter mixed convective boundary layer flow of 
nanofluid past a permeable stretching sheet in a Darcian porous 



media. The governing boundary layer equations are transformed 
to non-similar equations using suitable transformations. Numerical 
solutions are obtained using an implicit finite difference method. 
Validation with published paper is achieved. The effect of the 
controlling parameters on the dimensionless velocity, temperature, 
friction factor and heat transfer rates are illustrated via figures. 

Nanofluid Dynamic Transport Model 

Consider the steady 2-D laminar boundary layer flow of viscous 
water-based nanofluids containing Cu and AI2O3 nanoparticles. 
The base fluid and the nanoparticles are assumed to be in thermal 
equilibrium. It is assumed that the plate is moving with a velocity 
w^v = c X in the quiescent free stream. It is further assumed that the 
plate surface is subjected to linear hydrodynamic and thermal slip 
boundary conditions. We consider the nanofluid is a two 
component mixture with the following assumptions: (i) incom- 
pressible; (ii) no-chemical reaction; (iii) negligible viscous dissipa- 
tion; and (v) nano-solid-particles and the base fluid are in thermal 
equilibrium. The effect of g-jitter is induced by mixed convective 
flow of a nanofluid past a permeable vertical stretching. Following 
Sharidan et al. [46], the gravity acceleration is given by 

g*(/) = go[l +^ cos(n OJ t)] K , where go is the time-averaged value 
of the gravitational acceleration, g*(0 acting along the direction 

on the unit vector K, which is oriented in the upward direction, £ 
is a scaling parameter, which gives the magnitude of the gravity 
modulation relative to go? t is the time and co is the frequency of 
oscillation of the g-Jitter driven flow. If e < < 1 then the forcing 
may be seen as a perturbation of the mean gravity. The scheme of 
physical configuration is shown in Fig.l. Field variables are 
Darcian tangential and normal velocity components u, v and 
temperature T. Under these assumptions, the basic continuity, 
momentum and energy equations in dimensional form can be 
written as (Tiwari and Das [15], Vajravalu et al. [47]) 



dx dy 



V.V 

\, (i ) \bmentum boundaty la>er and 
(ii) Thennal boundarv* la>'er 




Figure 1. Geometry of the flow model and coordinate system. 

doi:10.1371/joumal.pone.0099384.g001 
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du du du 
— -\-u— +v— = 

ot ox oy 



1 d 

pnf^y 



du 



Pnf 



kpPnf 



U, 



(2) 



T = (jCft, n = A V, u = cx — , 

T-Tao , ... g*(/) 



(6) 



T —T 

w oc 



^0 



dT dT dT 



(3) 



The initial and boundary conditions are (Karniadakis et al. [32]) 
t<0, u=v = 0 any x, y 



Substituting Eqn. (6) into Eqns. (2)-(3), we obtain the following 
differential equations 



dff' df] df]^ df]^ \drj 

Da drj ^ dxdr] ' 



(7) 



t^O^ u = u,^-\-NiVnf ^,y= —Vw, T=T,^-\-Di ^ at y = 0, 
dy dy (4) 

u = 0, T^Tao as y^co, 

where p^y, (Xnf, knf are the properties of nanofluids and are 
given by (Oztop and Abu-Nada [48]) 



knf^ {ks + 2kf)-2^ (kf-ks) 
kf {ks + 2kf)+^ {kf-ks) ' 



{P^p)nf 



(5) 



Here py is the density, /Zy is the viscosity, kf is the thermal 
conductivity of the base fluid, (J) is the solid volume fraction 
parameter of the nanofluid = 0 corresponds to clear fluid), is 
the the density of the nanoaprticles, ks is the thermal conductivity 
of the nanoparticles, kp is the permeability of the porous media, 
N\ is the velocity slip factor having dimensions of second/meter, 
Di is the thermal slip factor having dimensions of meter, v^i; is the 
mass transfer velocity, Vvv > 0 is for suction and < 0 for injection. 
Further, {pCp)^^ is the heat capacitance of the nanofluid, Cp is the 
specific heat at constant pressure. Following Vajravalue et al. [47], 
the thermal expansion coefficient of the nanofluid can be 
determined by {pP)^f = {\—^){pP)f-\-^{pP)^. Following Mu- 
khopadhyay and Layek [49], Hamad et al. [44] we assume the 
temperature dependent viscosity varies according to 
/V= /V 00 [l+^o(^w-^)]=/V^ [1+^(1-^)], where ^f^ is 
the constant undisturbed viscosity is constant with /?o > 0, 0 is 
defined in Eqn. (6), A = b()(T^^, — Tco) is the viscosity variation 
parameter. 

To reduce the number of dependent variables as well as number 
of equations, we use stream function ij/ defined by 



dij/ 



V = — Note that xj/ satisfies equation of continuity 

automatically. Now, following Sharidan et al. (2006), we introduce 
following transformations 



(8) 



The boundary conditions (4) become 

5/ 



The constants Ci,C2,C3 are defined as: 

[1+^(1-0)1 



^1 



C2= [1 +£ C0S(7rT)], 



C3 = 



Kfjkf 

(\-^)^^[pCp)j (pCp)/ 



with 01 = (1-0) 



and 02 = 



(i-^)+0 



(1-0)+^- 



pf 



where Pr = is the Prandtl number, a = -JcVf^Nx is the 

velocity slip parameter, b = ^J~cJvf^D\ is the thermal slip 
CO . 

parameter, Q = — is the non-dimensional frequency, e is the 
amplitude of modulation, X = — 2 the mixed convection 



(t, 0) = 1 + aCi (t,0), /(t, 0) ^(t,0) = l+b- (t,0), 
(t,oo) = 0(t,oo) = 0. 



(10) 



(11) 



(12) 



(13) 



(14) 
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Table 1. Thermophysical properties of the base fluid and nanoparticles [48]. 





Physical properties 


Base fluid 


Nanoparticles 






Water 


Cu 


AI2O3 


Cp (J/kg-K) 


4179 


385 


765 


P(kg/m^) 


997.1 


8933 


3970 


k(W/m-K) 


0.613 


401 


40 


ax 10-^ (m^/s) 


1.47 


1163.1 


1738.6 



doi:10.1371/joumal.pone.0099384.t001 



parameter, Gfx = 2 local Grashof number, 



knC 



is the suction/ 



Da = — ^ is the Darcy number and = ^ 
v/00 V^cly- 

injection parameter, Re^ = — — is the local Reynolds number. 

Thermophysical properties of the base fluid and nanoparticles are 
reported in Table 1. 

For the first level of truncation (or steady state), the time 
derivatives in Eqns. (7)-(8) can be neglected. Thus the steady state 
governing equations for the first level of the truncation can be 
written as 



df]^ df]drf- drf- \drj 



Da drj 



Pr drj^ drj 



The boundary conditions (9) become 



(15) 



(16) 



^ (0) = 1 + aCi 0 (0), /(O) =/„., 0(0) = 1 + 6 £ (0), 



drj 
drj 



drj 



(00)= ^(oo) = 0. 



The steady state solution of Eqns. (15) and (16) with boundary 
conditions (17) is found numerically and compared with Grubka 
and Bobba [50] in Table 2. 

Assuming new functions F(T,rj),0(T,rj) which are defined by 



ox OT 



(18) 



and restoring all of the terms neglected in the first level of 
truncation, the governing equations can be written as 



drj^ drj drj^ drj^ \drj 

Da orj orj 



Pr drj'^ drj 



(19) 



(20) 



The boundary conditions (4) become 

5/ 



drj 

(17) drj 



(t, 0) = 1 + aCi ^ (t,0), /(t, 0) ^(t,0) = \+b — (t,0), 



drj'^ 

(t,oo) = ^(t,oo) = 0. 



drj ^ 



(21) 



Table 2. Comparison of dimensionless heat transfer rates — 


O'ifd) for different values of Pr when s = Q.=fw = a = b = (l) = A = X = 0. 




Pr Grubka and Bobba [50] 


Present Results 


0.72 0.4631 


0.46325 


1 0.5820 


0.58198 


3 1.1652 


1.16524 


10 2.3080 


2.30800 


100 7.7657 


7.76565 


doi:1 0.1 371/joumal.pone.0099384.t002 
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Cu-Water 




a=0, 1, 2 




Al203-Water 



0 - 



A 
0 
2 



b=Da=/.=T=0.5 
£= (})=0.1, Q=1 
f =0.3 



L(aL 



n 



_L(b)_ 



n 



Figure 2. Effects of viscosity variation and velocity slip parameters on dimensionless velocity for (a) Cu-water and (b) AlzOs-water 
nanofluids. 

doi:10.1371/journal.pone.0099384.g002 



Eqns. (19) and (20) contain two additional functions. In order to 
solve these equations, we need two more equations which can be 
obtained by differentiating Eqns. (19)-(20) with respect to t. 
Neglecting derivatives of F and 0 with respect to t, and 
simplifying, we get 



with new boundary conditions 



(23) 



Ci Ci 



0/ — ^;^77; F — C2 ASK sm {mjO ^ ^ 



DaC\ ^ DaC\ 

-C2A[1+£COS(71t)]0 = O 



Cu-Water 



0.6 1^=-0.3, 0, 0.3 



Da 

0.1 

1 



b=a=/-=T=0.5 
s= (j)=0.1, Q=A=1 




Aipj-Water 



0.6 



tt=-0.3, 0,0.3 




Da 
0.1 
1 



b=a=/=T=0.5 
s= (tF0.1,O=A=1 



0 '- ^ 



L(aL 



_l I I I I I I I I I I I I I I I L_ 

1 2 3 

Tl 



Figure 3. Effects of Darcy number and suction/injection parameters on dimensionless velocity for (a) Cu-water and (b) Al203-water 
nanofluids. 

doi:10.1371/joumal.pone.0099384.g003 
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Cu-Water 



0.6 



1^0.4 ^\ 



0.2 



L(aL 



({fO, 0.1,0.2 



1^. 



0.3 
-0.3 



'I 



A=Da=b=a=/.=i=0.5 
s=O=0 



_l I I I I I I I I I I I I I I I I I L. 



4 



Al203-Water 



(t)=0,0.1,0.2 




0.3 
-0.31 



A=Da=b=a=A=T=0.5 
£= Q=0 



_l I I I I I I I I I I I I I I I I I I 



4 

Tl 



Figure 4. Effects of suction/injection parameters and solid volume fraction of nanoparticles on dimensionless temperature for (a) 
Cu-water and (b) Al203-water nanofluids. 

doi:1 0.1 371 /journal. pone.0099384.g004 

F(0) = ^ {-A0(O)f "(0) + [1 + A(1-0(O)]F'(O)}, 

F(0) = 0,0(0) = Z70'(O),/"(oo) = O,0(oo) = O Pful' " kf{T„-T^)' ^ > 



where Tv^, qw are shear stress and the wall heat flux (dimensional) 
and are defined as 



Quantities of Physical Interest 



The quantities of enmneerinsf interest, in this study, are the local / du\ , / dT\ 



skin fi:'iction factor Cfx and the local Nusselt number Nu^ can be " ^"^ V^j/ ,=o ' ~ V ,=o 

found from the following definition 

Using Eqns. (6) and (26), we have from Eqn. (25) 
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Cu-Water 




Q 
0 

10 



Da=A=1,f^.=0 

/ =T=a=0.5, £=(})=0.1 



Al203-Water 




Q 

0 

10 



Da=A=1,f^=0 
/=i=a=0.5, £=(tF0.1 



-I I I I I I L_ 



J I I I l_ 



(b) 



n 



Figure 6. Effects of thermal slip and dimensionless frequency on dimensionless temperature for (a) Cu-water and (b) AlzOs-water 
nanofluids. 

doi:10.1371/journal.pone.0099384.g006 



kf or] 



where Re^^ = is the local Reynolds number. 

V/oo 



Validation with Published Results 

It is interesting to note that if we put 
A = a = b = X = (l) = 0, Da^co, in Eqns. (15)-(17), we have the 
same Eqns. as derived by Grubka and Bobba [50]. A comparison 

(27) 

of — 6'(0) is presented in Table 2 for the special case of steady- 
state flow (£ = Q = 0), no suction/injection (fw = ^), no slip 
{a = b = 0), conventional fluid (^ = 0) and constant viscosity 
(^ = 0). The present numerical results are found to be in good 
agreement for each value of Pr. We, therefore, are confident that 
our results are accurate. 



Cu-Water 




0.5 



L a=b=Da=T=1/2 4>=0,0.1,0.2 
/ =0, £=0=0.2 

_l I I I I I I I I I I I I I I I I I L_ 



-0.5 



0 



0.5 



2.5 



2 - 



0) . 



0.5 - 



Al203-Water 



A 

: 0 (t>=o, 0.1 , 0.2 






: 1 1 

- ^^^^ — .>--^ 


^ _^ ^ — * ^ 






: a=b=Da=i=1/2 




> =0, £=Q=0.2 

1 1 1 1 1 1 1 1 1 1 1 1 1 


1 1 1 1 1 1 1 



-0.5 



_L(">L 



0.5 



Figure 7. Variation of skin friction with suction/injection parameters and solid volume fraction of nanoparticles for different values 
of viscosity variation parameter for (a) Cu-water and (b) AlzOs-water nanofluids. 

doi:1 0.1 371 /journal. pone.0099384.g007 
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Cu-Water 



0.7 




M a=b=/=1/2, ({)=0.05 

- ^f/ Da=1,A=f^.=0 

I I I I I I I I I I I I 



0.5 



L(a)_ 



1 

T 



1.5 



0.98 



Al203-Water 



0.96 

X 



o 

CM 

^= X 

0) 
0^ 



0.94 - 



0.92 



Q 

0.1 

0.5 ^ 



s=0.2, 0.3, 0.4 




yffr a=b=/.=1/2, (t)=0.05 
/ Da=1,A=f,=0 

ftr I \^ Lj \—x \ L_i i_j i_l 



0.5 



'(b) 



1 

T 




Figure 8. Variation of skin friction with dimensionless time and scaling parameters for different values of dimensionless frequency 
for (a) Cu-water and (b) Al203-water nanofluids. 

doi:10.1371/journal.pone.0099384.g008 



Results and Discussions 

Equations (19), (20), (22) and (23) with boundary conditions (21) 
and (24) were solved numerically to obtain the dimensionless 
velocity and temperature fields as well as the skin friction and heat 
transfer rate. The effects of viscosity variation and velocity slip 
parameters on the dimensionless velocity are compared for Cu- 
water and Al203-water nanofluids in Figs. 2(a) and 2(b). The other 
parameters are kept fixed and suction is considered in this case. 
For uniform viscosity and no slip condition, the dimensionless 
surface velocity of both nanofluids is found to be higher in 
unsteady flow. But as slip parameter increases, the dimensionless 
velocity at the surface decreases and converges quickly. As a result. 



the hydrodynamic boundary layer thickness and hence the skin 
friction decreases. Physically, as the velocity slipping parameter a 
increases the differences between the wall and the fluid velocities 
near to the wall rises. Increasing the slipping factor may be looked 
at as a miscommunication between the source of motion (the plate) 
and the fluid domain. Note that the case a = 0.0 corresponds to 
conventional no slip boundary conditions. It is noticed from the 
same figure that the dimensionless velocity increases with an 
increase in viscosity parameter for both in the case of slip and no 
slip boundary conditions. Physically, with an increase in the 
viscosity parameter A, fluid viscosity decreases resulting the 
increment of velocity boundary layer thickness. Due to lower 
density of AI2O3 nanoparticles, the dimensionless velocity of 



2.4 
2.1 
1.8 



o 

0.9 



Cu-Water 



0.6 
0.3 



b=Q=0.5, 1=1 



Da 

1 

3 



ta=0, 0.5, 1 

4>=0.1, A=f^.=0,£=0.2 

I I I I I I I I I I I I I I I I I I I I I I I 



0.2 



0.4 0.6 



0.8 



Al203-Water 



2.1 
1.8 
1.5 



O 1.2 1 



0.9 



b=Q=0.5, 1=1 



Da 
• 1 
> 3 



0.6 

o.3|^a'=0, 0.5, 1 
0 



(t>=0.1,A=t=0,£=0.2 



I I I I I I I I I I I I I I I I I I I I I I I I 



0.2 



(b) 



0.4 0.6 



0.8 



Figure 9. Variation of skin friction with mixed convection and velocity slip parameters for different values of Darcy number for (a) 
Cu-water and (b) AlsOs-water nanofluids. 
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Al203-water nanofluids is found to be smaller at the surface. The 
Darcy number, in fact, shows the ability of fluids to flow through 
porous medium. The greater the Darcy number the greater will be 
the flow of fluids and hence the greater velocity of fluids. This is 
illustrated in Figs. 3(a) for Cu-water and in Fig. 3(b) for AI2O3- 
water nanofluids respectively. The suction velocity tends to reduce 
the boundary layer thickness whereas the boundary layer thickness 
increases in case of injection. This is shown in Figs. 3(a) and (b) for 
both nanofluids. Physically, suction causes the boundary layer to 
adhere more closely to the wall and this destroys momentum 
leading to a plummet in velocity. Hydrodynamic boundary layer 
thickness is therefore decreased with suction. Injection adds 
nanofluid velocity via lateral mass flux through the sheet and this 



assists momentum development, enhancing velocity and causing a 
concomitant increase in momentum boundary layer thickness. 
Again due to lower density of AI2O3 nanoparticles, the dimen- 
sionless surface velocity of Al203-water nanofluids is found to be 
lower. The effects of different controlling parameters on the 
dimensionless temperature are shown in Figs. 4-6 for both Cu- 
water and Al203-water nanofluids. The effects of solid volume 
fraction of nanoparticles on the dimensionless temperature are 
shown in Figs. 4(a) and (b) for suction and injection of both 
nanofluids. In unsteady flows, the solid volume fraction of 
nanoparticles increases the dimensionless temperature of both 
nanofluids inside the thermal boundary layer and as a result the 
thermal boundary layer thickness increases for the selected values 
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of the controlling parameters. The suction/injection parameters 
have strong effects on the thermal boundary layer thickness. In 
case of suction, the thermal boundary layer thickness is found to be 
smaller than in injection for both nanofluids. However, due to 
higher thermal conductivity of Cu-water nanofluids, the dimen- 
sionless surface temperature of Cu-water nanofluids is found to be 
lower than Al203-water nanofluids. The effects of the dimension- 
less time and the scaling parameter on the dimensionless 
temperature are depicted in Figs. 5(a) and (b) for both nanofluids. 
The scaling parameter e shows the magnitude of the gravity 
modulation relative to go- It can be seen that for steady case, when 
T = 0, the surface temperature is found to be lower for smaller 
values of e in both cases. As the dimensionless time increases, the 
surface temperature increases for both nanofluids. However, both 
the dimensionless time and the scaling parameter have negligible 
effect on the thermal boundary layer thickness. The variation of 
the dimensionless temperature with the transverse distance inside 
the thermal boundary layer is shown in Figs. 6(a) and (b) for both 
nanofluids. The effects of thermal slip and dimensionless frequency 
are investigated while all the other parameters are kept constant in 
both cases. It can be seen that, in the absence of thermal slip, the 
surface temperature is highest and it reduces with an increase in 
the thermal slip. In steady-state (Q = 0), the dimensionless 
temperature is found to be higher in the thermal boundary layer 
in both cases. The variation of skin friction with different 
controlling parameters is shown in Figs. 7-9 for the selected 
nanofluids. The skin friction is found to increase in case of suction 
and decrease in case of injection as shown in Fig. 7(a) for Cu-water 
and in Fig. 7(b) for Al203-water nanofluids respectively. The 
effects of temperature dependent viscosity on skin friction can also 
be observed in these figures. When the viscosity is constant, the 
skin friction is found to be lower and it increases with temperature. 
It is important to note that as we go from suction to injection 
region, the effects of viscosity variation parameter diminish. In the 
absence of nanoparticles (pure water), the skin friction is found to 
be lowest. But as the solid volume fraction of nanoparticles 
increases, the skin friction increases due to higher density of 
nanofluids. Since the density of Cu-water nanofluids is higher than 
Al203-water nanofluids, the skin friction for Cu-water nanofluids 



is found to be higher. Due to sinusoidal nature, the skin friction 
first increases up to maximum value and then decreases with 
dimensionless time. This is shown in Figs. 8(a) and (b) for two 
selected water-based nanofluids when there is no suction or 
injection. It can be seen that the scaling parameter has no 
appreciable effect on the skin friction in any case. The variation of 
skin friction with mixed convection parameter X for different 
values of Darcy number and momentum slip is shown in Figs. 9(a) 
and (b) for two selected nanofluids. In case of forced convection 
(/l = 0), the skin friction is found to be higher than in free 
convection (A=l) in both case. Both Darcy number and 
momentum slip tend to reduce skin friction in both cases. 
Figures 10-12 show the variation of Nusselt numbers with 
different controlling parameters for two selected water-based 
nanofluids. It can be seen that the Nusselt numbers for pure water 
((/> = 0) are lower and increase with an increase in the solid volume 
fraction of nanoparticles, as shown in Figs. 10 (a) and (b) for both 
nanofluids. This is due to increase in thermal conductivity with the 
solid volume fraction of nanoparticles. However, the effects of 
viscosity variation parameter on the Nusselt numbers are not 
appreciable. Like skin friction, Nusselt numbers also increase 
during suction and decrease during injection. The sinusoidal 
nature of Nusselt numbers can be observed in Figs. 1 1 (a) and (b) 
for both nanofluids. The case of uniform viscosity without any 
suction/injection is considered for investigating the effects of 
scaling parameter and dimensionless frequency on the Nusselt 
numbers for both nanofluids. The effects of mixed convection on 
the Nusselt numbers for different values of Darcy number and 
momentum slip are shown in Figs. 12(a) and (b) for both 
nanofluids. As expected, the Nusselt numbers increase with mixed 
convection parameter and Darcy numbers. But this increase is 
found to be negligible in no slip condition. As the momentum slip 
increases, Nusselt numbers decrease. Since the thermal conduc- 
tivity of Cu is higher than AI2O3, the Nusselt numbers are found to 
be higher for Cu-water nanofluids. 
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Conclusion 

The two-dimensional unsteady laminar g-jitter mixed convec- 
tive boundary layer flow of Cu- water and Al203-water nanofluids 
past a permeable stretching sheet in a Darcian porous is studied by 
combined non-similar and numerical solution technique. Main 
findings are given below. 

(1) Dimensionless velocity reduces as velocity slip parameter 
increases for both nanofluids and for fluid with constant and 
variable properties. 

(2) Dimensionless velocity increases as Darcy numbers increase 
for both nanofluids in the case of both permeable and 
impermeable plate. 

(3) Dimensionless surface temperature increases as nanoparticle 
volume fraction increases for both nanofluids 
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